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We have introduced a mechanism for synchronizing chaotic systems by one way couphng with a 
variable delay that is reset at finite intervals. Here we extend this method to time delay systems and 
suggest a new cryptosystem based on this. We present the stability analysis as applied to time delay 
systems and supplement this by numerical simulations in a standard time delay system like Mackey 
Glass system. We extend the theory to multi- delay systems and propose a bi channel scheme for 
the implementation of the scheme for communication with enhanced security. We show that since 
the synchronizing channel carries information from transmitter only at intervals of reset time, it 
is not susceptible to reconstruction. The message channel being separate can be made complex 
by linear combination of transmitter variable at different delay times using mutiple delay systems. 
This method has the additional advantage that it can be adjusted to be delay or anticipatory in 
synchronization and these provide two additional basic keys that are independent of system delay. 

PACS numbers: 05.45.Xt,05.45.Vx 



I. INTRODUCTION 

In recent years, a number of cryptosystenis based on 
synchronized chaotic systems have been introduced [l|,l2|- 
Many of these systems generally make use of the broad 
band spectrum of the chaotic output from the trans- 
mitter as a carrier for secure masking of the signal and 
its synchronization with the receiver for easy recovery 
of the signal. In addition to chaotic masking, several 
other schemes involving chaotic modulation and chaotic 
switching etc. are also used for a practical realization of 
the chaos based cryptosystem 0-|3- It has been shown 
that message transmission using synchronization can be 
implemented in practice using semi conductor lasers and 
electronic circuits of chaotic systems piil • The most fre- 
quently studied schemes involve complete synchroniza- 
tion of one way coupled identical systems. However 
generalized synchronization of non identical systems and 
phase synchronization have been adapted to secure com- 
munication using Chua's circuits [a, |lO|. 

Of late, it has been shown that low dimensional chaotic 
systems can be easily unmasked by reconstructing the dy- 
namics from time series available from the channel by an 
intruder. Since then, the use of hyper chaotic systems 
derived from time delay systems have been proposed for 
improved security |11| . Their methods of implementation 
using time delayed optical systems and circuits and gen- 
eration of secret keys have also been studied in dctail[12- 
|2]| . However, later studies have indicated ways of recon- 
struction of such systems from time series and extraction 



'Electronic address: 
T Electronic address: 



g.ambika@iiserpune.ac.m| 



aniritkar@prl.res.in| 



of messages masked by hyper chaotic carriers |22h27I | . In 
general, such cracking schemes, are aimed at identify- 
ing the time delay parameter from the cipher text in the 
channel. For enhanced security, several variations like 
multiple delays, variation and modulation of de lays, ha r- 
monics and multi channels have been reported |28l43l| . 

We have recently introduced a mechanism for synchro- 
nization of chaotic systems by coupling involving a vari- 
able delay that is reset at finite intervals [32|. In this 
method, the delay varies in the same way as the system 
in time and hence the value of the synchronization signal 
remains constant for the reset interval at the end of which 
it is reset to the value at that time. We have applied 
this to standard low dimensional systems like Rossler and 
Lorenz and carried out the stability theory and Lyapunov 
exponents analysis for such a coupling scheme. 

In this paper, we further develop this method to apply 
to hyper chaotic time delay systems like Mackey Glass 
and Ikeda systems with possible applications to secure 
communication. In this context, we note that several 
studies on variable and multiple time delays in synchro- 
nized time delay systems have been reported |33| with ap- 
plications to communication |34l437| . But in our scheme, 
unlike the above methods, the variability is in the de- 
lay or anticipation time in the coupling function and not 
in the feedback delay in the system dynamics. We fur- 
ther suggest a bi-channel scheme for its implementation 
where the synchronization channel carries information 
about the transmitting dynamics only at intervals of re- 
set time and hence is not susceptible to reconstruction by 
an intruder. The message channel being separate can be 
made complex enough. At large reset times, the reported 
methods of reconstructing the dynamics by first comput- 
ing the delay time of the system dynamics is not possible, 
since the information of the system variable is available 



only at intervals of the reset time. By a suitable choice 
of the carrier as a combination with multiple delays from 
the transmitter for effective masking, we show how the 
message can be retrieved truthfully at the receiver once 
it synchronizes with the transmitter. 

The method reported here has the additional advan- 
tage that it can be easily adjusted to be delay, antici- 
patory or isochronous in synchronization. Anticipation 
in the context of communication is not essential even 
though it will give more time for decoding at the re- 
ceiver. However, the delay and anticipatory times serve 
as additional keys. In this paper we give more emphasis 
on small anticipatory times so that the reset time can be 
maximized with large enough system delay. This makes 
the method highly cost effective with more complex out- 
put for enhanced security. Such optimization is possible 
due to the ffexibility among the four different time scales 
in our method viz. delay in the transmitting system , 
anticipation in the receiving system, system delay and 
reset times. We present the stability analysis as applied 
to time delay systems and supplement this numerically 
by computing the stability regions for synchronization 
in the parameter space for a standard time delay sys- 
tem normally used for communication viz. Mackey Glass 
(MG) system. We also present the study of two MG 
systems with multiple delays using the present scheme. 
The applicability for communication is illustrated by a 
bi-channel scheme for an analogue signal. We establish 
the enhanced security of this method by explicitly apply- 
ing the existing methods of cracking applicable to such 
time delay systems. 



II. SYNCHRONIZATION OF TIME DELAY 

SYSTEMS WITH VARIABLE DELAY AND 

RESET 

We consider systems whose intrinsic dynamics follows 
a time delay equation 



the following dynamics for the coupled systems under 
this scheme. 
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where Xt-^ = x{t — ti), yt^ — y{t — t2), t^ is the reset- 
ting time and X(t'.t") is an indicator function such that 
X(t'.t") — I ior t' < t < t" and zero otherwise. Both 
the delays ii and f2 depend on time and we choose this 
dependence as 



ti 



Ti + t — rriT, 1 = 1,2. 



This takes care of the variation in the delay times men- 
tioned above. The synchronization manifold for the 
coupled systems ^ is defined hy y{t — T2) = x{t — 
Ti) or y{t) = x{t -Ti+ T2) 

We have shown that under this type of coupling, it 
is possible to have isochronal, delay or anticipatory syn- 
chronization between chaotic systems j32l |. In all these 
cases, the transverse system can be defined by A = 
y — Xri~T2- Its dynamics in linear approximation can 
be derived from Eq. ^ as 



A = -kA + f {Xri-rz+rJ'^^r, 
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where Am — A(i — ^2) = A{'niT — T2) is a constant in 
any reset interval. ti can eliminated by shifting the time 
scale of drive system and redefining T2 . 

We approximate the equation by replacing / by an 
effective time average value ly 



/(a^r,; 



(1) 



where Xr^ = x(t — Tg); t^ being the time delay in the 
system. 

Two such systems are coupled in the drive response 
mode with a linear difference coupling where the drive 
variable is delayed by ti and the response by T2. These 
delays vary in time with the same time scale as the sys- 
tem within an interval called the reset time r^. At the 
end of this interval, the delay is reset to the initial value 
and the dynamical variable set to its value at that time. 
This makes the coupling function constant over r^ and 
hence the information from the transmitter for the syn- 
chronizing signal is required only at intervals of r^ . The 
dynamics of the two coupled systems thus evolves under 
four different time scales in addition to the system dy- 
namics, viz the delay time ti, the anticipatory time T2, 
the reset time Tr , the system delay time Ts ■ We prescribe 



A = -kA -I- I^Ar^ ^ <^ X! X(mr,,.(m+l)r,) Am (4) 

m=0 

Assuming a solution A = e'^* for the equation, 

A = -kA + l^Ar^ (5) 

we get 

X= -K + ve-^^^ (6) 

Then the solution of Eq. ^ in any reset interval , 
TOT <t< {m + 1)t, is 

A = aAm + C„e^* (7) 

where a = e/A is the normalized dimensionless coupling 
constant, and Cm is an integration constant, which can be 
determined by matching solutions across the reset points. 




FIG. 1; The variation in A as the system delay increases for 
two cases, K < //(upper curve) k, > //(lower curve). 



Then the stability analysis detailed in [S^l , can be easily 
adapted to time delay systems also, the only difference 
being that here A depends on the system time delay Ts 
and has to be evaluated for each case using Eq. ^ . This 
transcendental equation can be solved to show that for 
J/ > and K < V, \ is positive and as Ts increases A de- 
creases, while for k > //, A is negative and as Tg increases 
A increases. This is shown in Fig. [TJ 

In this paper, since we are interested in applications to 
communication, T2 need not be large and we concentrate 
in the regime < r2 < t^. Then the upper limit of 
stability in the parameter plane T2/T — a as derived in 
[32| has a peak at T2p. For T2 < T2p, the upper limit of 
stability is given by 
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(8) 



while for larger values of T2 (T2p < T2 < t^.) it is given by 

"« = 1 _ e-Ar. (9) 

Moreover, in the present context, the maximum reset 
time for a given T2 may be more relevant and from our 
earlier analysis, it is 
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A similar analysis will furnish the stability criteria with 
A replaced as 



N 
A = —K + 2^ CLiVie '' 
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(12) 



The above analysis is supported by numerical simula- 
tions using a standard time delay system, Mackey Glass 
system, in the next section. 



III. MACKEY GLASS SYSTEMS WITH 
VARIABLE DELAY IN COUPLING 

The Mackey Glass system, first introduced as a model 
for blood generation in patients with leukemia, is well 
studied as a model exhibiting hyper chaos [3^ ■ The syn- 
chronization in two such coupled systems has been re- 
ported for various types of synchronization [39|. It has 
been applied extensively in the context of communica- 
tions, especially with its circuit equivalents [40l|. The 
synchronization regimes of two such systems with multi- 
ple delays have been reported recently [41[. So also cou- 
pling with multiple delays [42, E^ has been proposed for 
more secure communications. In this section, we apply 
the scheme of variable delay described in the last section 
to two coupled Mackey Glass systems. The dynamics of 
a single system is 



ATy 
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for given T2 and for 



Using Eq. dH) and Eq. dTU]), Trmax 
different Tg can be calculated. 

We further extent this scheme to systems with multi- 
ple delays that are more realistic models in many appli- 
cations. They occur in many situations, an example for 
such multi-delay systems being the semiconductor laser 
with multiple opto electronic feed backs, where additional 
feed back is often required to stabilize the output. More- 
over such systems are of practical importance especially 
in the context of communication where the higher com- 
plexity can lead to better message security. The dynam- 
ics for two such systems coupled using the present scheme 
is 



l+x{t-TsY 



(13) 



As reported in (38|, for a=l, b=2 and c==10, as the 
delay time r^ is varied, the system has a fixed point at- 
tractor for Tg < 0.471, a limit cycle for 0.417 < r, < 1.33, 
a period doubling sequence in the range r^ = [1.33, 1.68] 
and chaos for Tg > 1.68. Thereafter the number of posi- 
tive Lyapunov exponents increases linearly with Tg, and 
the system shows hyper chaotic behavior. In our simu- 
lations here, we choose a similar set of parameter values 
with Tg between 2.5 and 5.0 corresponding to the hyper 
chaotic region and couple two such systems using our 
scheme given in Eq. [5] as 
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FIG. 2: The nature of variation of the coupling function hav- 
ing variable delay with reset (full line) and the synchronization 
error(dashed curve) for two coupled Mackey Glass systems. 
Here Ts^2.5, r^^ 0.5, ri = .02, T2= 0.1 and e = 2.0 



FIG. 3: The limits of stability of the synchronized state for 
two coupled Mackey Glass systems in the parameter plane 
T2 — e. The dotted and dashed curves for the upper boundary 
correspond to the curves obtained by fitting Eq. [5] and Eq. |5] 
obtained from the stability analysis, while the points corre- 
spond to the values obtained by direct numerical simulations 
using the correlation index. Here Ts=2.5, Tr= 0.2 and ri = .02. 
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For a typical choice of parameters, r^ = 2.5, Tj.= 0.5, 
Ti=.02, T2= 0.1, and coupling strength, e = 2.0, we 
present the variation in the coupling function accord- 
ing to our scheme and the corresponding error function 
in Fig. [21 It is clear that the coupling function is not 
changing within each reset interval and the synchroniza- 
tion error approaches zero. 

By numerically integrating the coupled systems in 
Eq. [T31 we calculate the correlation coe fficient C =< 
yi{t)xi{t + T2) > /y/< xfit) X yl{t) > between xi(t) 
and yi{t) shifted by the effective T2 = \t2 — ti\ (hereafter 
referred to as T2 itself). The region of stability of the syn- 
chronized state is isolated as the region where C — 0.999 
and boundaries of stability fixed when C goes below this 
value. The stability region thus obtained in the parameter 
plane T2 — e is plotted in Fig. |3l The region below the up- 
per curve and above the lower curve corresponds to total 
synchronization with correlation function 1.0 . The dot- 
ted and dashed curves for the upper boundary correspond 
to the curves obtained by fitting Eq.|8]and Eq.[9]obtained 
from the stability analysis, while the points correspond 
to the values obtained by direct numerical simulations 
using the correlation index. The agreement is found to 
be very good. 

The maximum reset time for a given r^ obtained by 
numerical calculations is given in Fig. |4]for two different 
values of Tg, 2.5 and 5.0. It is interesting that as r^ in- 
creases, maximum reset time also increases for a given 
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FIG. 4: The maximum reset time for stability for a given Ts 
for two coupled Mackey Glass systems as a function of T2 for 
two values of Ts, 2.5 and 5.0. 



T2. This is a desirable feature in applications like com- 
munication. 

We extend our numerical computations to MG systems 
with multiple delays. The dynamics of two such systems 
with four delays is given by 
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The stability region in parameter plane, T2 — e obtained 
by direct numerical analysis of Eq. I15[ using the same 
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FIG. 5: The stability boundary for two coupled Mackey Glass 
systems with 4 delays in the parameter plane T2 — e obtained 
by direct numerical simulations using the correlation index. 
Here r,.= 0.5 and ri = .02. 



correlation index as above is given in Fig. [5] Here (6^, Tg.) 
values used are (2.0,2.5), (0.2,2.2), (0.05,2.3)and (0.1,2.4). 
We note that the region of stability in this case is more 
compared to single delay which may be useful in many 
applications. 



IV. APPLICATION TO SECURE 
COMMUNICATION 

In this section we propose a scheme for secure commu- 
nication using hyper chaotic systems using the method 
of synchronization discussed in this paper. We also show 
how the method can lead to enhanced security. We follow 
the chaotic masking method as applied to a bi-channel 
scheme. The transmitter and receiver dynamics in our 
scheme are given by Eq. [131 Here we take {bi,Ts^) as 
(2.0,5.0), (0.2,4.0), (0.05,3.0), (0.1,2.0). The synchroniz- 
ing signal S(t) is the intermittent values of x(t- ti) sent 
at intervals of the reset time r^ . 

The message to be transmitted is assumed to be a sine 
signal given by 



m{t) = 0.01sm(a;i) 



(16) 



This message is encrypted using method given in |41| 
for bi-channel scheme. The carrier is generated by com- 
bining the output of the transmitter at eleven different 
delay times using a nonlinear function as given below. 
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Here the values of {ai,Ts-) used are (4.0,2.0,), (5.0,1.9), 
(6.0,1.8), (7.0,1.7), (8.0,1.6), (9.0,1.5), (10.0,1.4), 
(11.0, 1.3), (12.0, 1.2), (13.0,1.1), (14.0,1.0). This carrier 



FIG. 6: The cipher text CT(t) given in Eq[T8] generated by 
masking the signal m(t) using carrier wave given in Eq. 1171 
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FIG. 7: The plain text PT(t) or recovered signal using Ea ll9l 



is used to mask the message and the cipher text C(t) 
thus obtained is 



C(t) = m{t) 
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(18) 



This cipher text is communicated to the receiver along 
one channel, while the synchronizing signal S(t) is trans- 
mitted through another channel. Once the systems at 
both ends are synchronized, the signal is recovered at 
the receiving end as plain text. 
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We display the cipher text C(t), the recovered plain text 
P(t) and the error in the recovered plain text in Fig. |6l 
Fig. [7] and Fig. [8] computed using the method. It is seen 
that after the transients, the error decreases to zero and 
the message can be recovered truthfully. 

To illustrate the enhanced security of the method, we 
note that S(t), derived using the present scheme, has in- 
formation about the dynamics of the transmitter only at 
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FIG. 8: The error in the recovered plain text as compared to 
the message m(t). The inset shows the Log\error\ for longer 
times. 
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FIG. 9: The method reconstruction in 1231] applied to the 
synchronizing signal to extract the system delay Ts- It is clear 
that there is no single absolute minimum. 



intervals of rest time r^ and as such can not be used to de- 
rive the dynamics using conventional methods of cracking 
reported so far. We show this explicitly by considering 
the method used in [22, [23] , where the time delay in the 
dynamics of the systems is to be extracted first to re- 
construct the transmitter system. For this the number 
of pairs of extrema separated in time by different r in 
the time series obtained from the channel is taken and 
the absolute minimum gives a clue to the time delay in 
the system t^ . Once this is obtained the functions in the 
transmitter dynamics can be derived [2J|. For applying 
this to our method, we reconstruct the time series from 
the synchronizing channel by removing the reset intervals 
and use this collapsed time series to take the count of ex- 
trema for different separations. This is given in Fig. [HI It 
is clear that as the minima happen at a large number of 
points due to finite reset times, an estimate of the time 
delay Ts cannot be obtained. 

We also use the same technique for checking possibility 
of reconstruction from the cipher text. Here for the same 



FIG. 10: The method of reconstruction in f2s| applied to the 
cipher text. The complexity of the cipher text gives several 
minima and hence it is dificult to get an indication of Ts . 



the cipher text in Eq. I18[ we find that the information 
about the system delay is sufficiently hidden. (Fig. [6|). 
Moreover, the eleven different sets of parameters (5i, Tg.) 
used in constructing the cipher text gives added secu- 
rity, by providing extra keys. In addition to this in our 
method, the two time delays ti and T2 serve as additional 
keys. 



V. CONCLUSION 

We have introduced a coupling scheme for synchro- 
nizing two one way coupled chaotic systems where the 
delay in coupling varies with system dynamics within in- 
tervals of the reset time. In this paper, we have extended 
this scheme to include time delay systems where the syn- 
chronization can be delay or anticipatory independent of 
the system delay. The stability analysis is adapted for 
such systems and applied to a standard time delay sys- 
tem like Mackey Glass system. We suggest the possibility 
of applying this method in communication with enhanced 
security. For this, we propose a bi channel method for 
implementation of the same. Our method has the advan- 
tage that synchronizing channel carries only intermittent 
information from the driver. We show that this enhances 
the security of transmission as extracting system delay 
from this channel is not possible by the conventional 
methods reported so far. The cipher text channel can 
be made complex resulting in extra keys. We illustrate 
this by using a carrier wave constructed from a multi de- 
lay Mackey Glass system as driver with different delays. 
Such bi channel schemes have been reported earlier also. 
But our scheme promises extra security in the synchro- 
nizing channel due to intermittent transmission and the 
delay and anticipation in synchronization in our method 
function as additional keys that are independent of sys- 
tem delay. 
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